Preprint typeset in JHEP style - HYPER VERSION 



December 2012 



Positroids, Plabic Graphs, & 
Scattering Amplitudes in Mathematica 



Jacob L. Bourjaily 

Department of Physics, Harvard University, Cambridge, MA 02138 

Abstract: The many intricate connections between scattering amplitudes, on-shell 
diagrams, and the positroid stratification of the Grassmannian has recently been de- 
scribed in detail. In order to facilitate the exploration of this rich correspondence, we 
have prepared a public Mathematica package called "positroids" which includes 
an array of useful tools including those for the construction of canonical coordinates 
for positroid configurations, the drawing of representative on-shell (plabic) graphs, 
and the evaluation of on-shell differential forms. This note documents the functions 
made available by the positroids package; the package's source code together with 
a Mathematica notebook containing many detailed examples of its functionality 
are included with this note's submission files on the arXiv. 
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1. Introduction 

The recent work of [1] presents a comprehensive summary of the extensive cor- 
respondence between "on-shell diagrams", [2-6], the Grassmannian contour inte- 
gral described in [7-11], scattering amplitudes in planar, maximally supersymmetric 
("A/" = 4") Yang-Mills (SYM), [12-18], and the combinatorics and geometry of what 
is known as the positroid stratification of the Grassmannian, [19,20]. At the heart of 
this story is the fact that scattering amplitudes can be represented (to all loop orders) 
in terms of on-shell diagrams, and that (reduced) on-shell diagrams can be fully char- 
acterized combinatorially by permutations — associated with left-right paths; more- 
over, these same permutations label the configurations of the positroid stratification 
of the Grassmannian G(k, n) of /c-planes in n dimensions. These strata are naturally 
endowed with positive coordinates on and a canonical volume-form, [21,22], which 
when expressed in terms of positive coordinates, is simply: rflog(«i) A- ■ ■Adlog(ad). 
Because on-shell diagrams can be directly represented (and computed) as integrals 
of this invariant volume-form over their corresponding positroid configurations, this 
makes the evaluation of on-shell diagrams exceedingly simple. 
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We can illustrate this rich correspondence with the following example: 




1 a$ (a 5 + a 14 a 8 ) a 5 au 











I) 







1 aio (ai ai3+Q4) 0407 

1 Q'g (a3 + "l2Q'6) j 

aian — ctict2 — aiQf2d7 1 



e G(4, 9) 



t 



h 2345678 9\ 

14* 4- 4- 4- 4- 4- 4- 4- 4-1 

\7 5 6 1 8 9 4 2 3/ 

a={7, 5, 6, 10, 8, 9, 13, 11, 12} 



t 



3 



2 




da 



A- • -A 



14 S kx4 (C(a)-?j)S kx2 (C(a)-\)6 2x{n - k) (X-C(a) ± ) 



Starting with the on-shell diagram on the left, we find that it would be labeled 
(via left-right paths) by the permutation denoted a = {7, 5, 6, 10, 8, 9, 13, 11, 12}; this 
permutation also labels the Grassmannian configuration drawn (projectively 1 ) on the 
right — a configuration represented by the matrix C(a) above, parameterized by the 
positive coordinates in terms of C(a), the corresponding on-shell 'function' Jo- 
associated with the diagram is determined by the integral at the bottom. 

We will not review these ideas here, but instead refer the interested reader to 
reference [1] for a thorough introduction and summary together with a more com- 
prehensive list of references to the existing literature. 

In order to help facilitate further investigation along these lines, however, we have 
prepared a Mathematica package called "positroids" — which is documented in 
this note. The positroids package makes available many of the essential tools 
required to investigate the myriad connections between on-shell physics, scattering 
amplitudes, and the combinatorics and geometry of the positroid stratification of the 
Grassmannian. 



^^Hcre, each dot represents a column of the matrix C{a) viewed projectively as a point in P 3 . 



2. The MATHEMATICA Package positroids 



2.1 Obtaining the positroids Package and Demonstration Notebook 

From the abstract page for this paper on the arXiv, follow the link "other formats" 
(below the option for "PDF") and download the "source" for the submission. The 
source-file will contain 2 the positroids package's main source-code (positroids .m), 
together with a demonstration notebook (positroids_package_demo .nb) which de- 
scribes with detailed examples many of the functions defined by the package. 

2.2 Using the positroids Package 

Upon obtaining the package, users should open and evaluate the Mathematica 
notebook postroids_package_demo .nb; this notebook will copy positroids. m to 
the user's Application directory for MATHEMATICA; this will allow positroids to 
be started in any future notebook via the command: 

in[i]:= <<positroids .m 



(If the file "positroids .m" has not yet beed copied to the user's Application direc- 
tory, then the package can be initialized by saving a notebook to the directory where 
positroids .m is located, and evaluating "SetDirectory [NotebookDirectory [] ] " 
prior to the command "<<positroids .m".) 

2 Occasionally, the "source" file downloaded from the arXiv is saved without any extension; this 
can be ameliorated by manually appending ".tar.gz" to the name of the downloaded file. 



Out [1] 
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3. Functions Defined by positroids 

3.1 Operations on Permutations Labeling Positroids 

• boundary[permutation_ ]: returns a list of permutation labels for positroid cells 

in the boundary, d, of the cell labeled by permutation. For example, the co- 
dimension one boundaries of the generic configuration in G + (3,6) are: 

in[i]:= boundary [{4,5,6,7,8,9}] 

Out[i]:= {{5, 4, 6, 7, 8, 9}, {4, 6, 5, 7, 8, 9}, {4, 5, 7, 6, 8, 9}, 
{4, 5, 6, 8, 7, 9}, {4, 5, 6, 7, 9, 8}, {3, 5, 6, 7, 8, 10}} 

• che ckOper at or \_permutation_ ] \_operator_ ] : returns: 

• if the minor (operator) vanishes for the cell labeled by permutation; 

• 1 if minor (operator) is non-vanishing for the cell labeled by permutation. 

in[i]:= {checkOperator [{3,5,6,7,8, 10}] [{1,2,3}] , 
checkOperator [{3,5,6,7,8,10}] [{2,3,4}] } 
0ut[l]:= {0,1} 

• cyclicize[permutation_] : returns a sorted list of non-repeating permutations in 

the same cyclic class as permutation. 

in[i]:= cyclicize [{6,5,8,7, 10,9, 12, 11}] 

outci]:- {{4, 7, 6, 9, 8, 11, 10, 13}, {6, 5, 8, 7, 10, 9, 12, 11}} 

• cyclicRep[permutation_] : returns the lexicographically-first permutation in the 

cyclic-class of permutation. 

• decorate[permutation_] : takes an 'ordinary' permutation a of n integers, and re- 

turns a decorated, affine permutation a, adding n to the image of any a such 
that a(a) < a. For example, applying decorate to the ordinary permutation 
associated with the on-shell diagram given in the introduction (section 1) gives: 

in[i] : = decorate [{7,5,6,1,8,9,4,2,3}] 
out[i]:= {7,5,6,10,8,9,13,11,12} 

• dimension[permutation_] : gives the dimension of the positroid stratum labeled by 

permutation; for example, 

in[i]:= dimension [{3, 5, 7, 6, 8, 14, 10, 12, 11, 13, 16, 21}] 

0ut[l]:= 20 
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• dualGrassmannian[permutation_ ]: if permutation labels a (2?7,-4)-dimensional cell 

in the 'momentum-space' Grassmannian C &G(k+2,n), then dualGrassmannian 
returns the permutation label of the corresponding, 4/c-dimensional cell in the 
'momentum- twistor' Grassmannian C &G(k,n), and vice-versa; e.g., 

in[i] := dualGrassmannian [{6, 5, 8, 7, 10,9, 12, 11}] 

dualGrassmannian@% 
out[i]:= {2,5,4,7,6,9,8,11} 

{6,5,8,7,10,9,12,11} 

• eulerCharacteristicTable[n_,if_] : returns a table indicating the numbers of d- 

dimensional cells in the positroid stratification of G + (k,n). (This data is gener- 
ated using the combinatorial results of [20].) 

in[i] := eulerCharacteristicTable [10 , 5] 



dim 


# 


# 


dim 


24 


1 


10 


25 


22 


55 


220 


23 


20 


715 


2002 


21 


18 


4985 


11240 


19 


16 


23210 


44220 


17 


14 


78087 


128100 


15 


12 


195315 


276450 


13 


10 


362175 


437112 


11 


8 


482670 


482940 


9 


6 


432060 


339360 


7 


4 


228102 


126420 


5 


2 


54600 


16800 


3 





3150 


252 


1 




1865126 
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Euler Characteristic: 1865126 - 1865125 = 1 



• intersectionNumber[permutation_,m_:4 ] : returns the number of isolated solu- 
tions to j2]-dimensional kinematical constraints, T m (C), where C is the positroid 
labeled by permutation. Supposing that C is an (m fc)-dimensional cell in the 
momentum-twistor Grassmannian 3 , intersectionNumber counts the number of 

3 interscctionNumbcr also tests (2n— 4)-dimensional cells in the momentum-space Grassmannian. 
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isolated points C* EC f] Z 1 , where Z is a generic configuration in G(m, 



n 



In[l] : = 



Out[l] : = 



{intersectionNumber[{6,5,8,7, 10,9, 12, 11}] , 
intersect ionNumber [{15 , 14 , 8 , 7 , 21 , 20 , 19 , 13 , 12 , 26 , 

25,24,18,17,31,30,29,23,22,36}]} 

{2,34} 



■ inverseBoundary[permutation_] : returns a list of permutation labels for positroid 

cells in the invers e-boundary, d~ l , of the cell labeled by permutation. 

< legalPermQ[permutation_]: tests whether the list permutation does in fact denote 
a (possibly decorated) permutation on n integers. 

■ necklace[permutation_ ] : if permutation labels a positroid configuration C G G(k, n) 

then necklace returns a list of n, fc-tuples = [A\ , . . . ,A^) denoting the 
lexicographically- minimal non-vanishing minors starting from each column a. 
Recall from [1,19] how the necklace encodes the list of all ranks of consecutive 
chains of columns of C; for example, the necklace for the configuration in G(4, 8) 
labeled by the permutation {3, 7, 6, 10, 9, 8, 13, 12} would be given by, 

•13- 




AW 


= (8 91013) 


A (7) 


= (78910) 


AW 


= (6 7910) 




= (5 6 710) 


A W 


= (4567) 


A ® 


= (3457) 


A m 


= (2345) 


A m 


= (1245) 



(3.1) 



This data is generated by the function necklace according to: 



in[i] := necklace [{3,7,6,10,9,8,13, 12}] 
0ut[i]:= {{1,2, 4, 5}, {2, 3, 4, 5}, {3, 4, 5, 7}, {4, 5, 6, 7}, 
{5, 6, 7, 2}, {6, 7, 1,2}, {7, 8, 1,2}, {8, 1,2, 5}} 
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necklaceR[permutation_]:if permutation labels a positroid configuration C G 
G(k,n), then necklace returns a list of n, fc-tuples A^ a ^ = , . . . , A%) de- 
noting the lexicographically- maximal non-vanishing minors starting from each 
column a. Like the more familiar ' Grassmannian necklace' generated by the 
function necklace, this data similarly encodes the ranks of all consecutive chains 
of columns of C; to see this, consider the 'reverse' necklace for the configuration 
in G(4, 8) labeled by the permutation {3, 7, 6, 10, 9, 8, 13, 12}, 






= (4567) 


1(7) 


= (2 4 5 6) 


1(6) 


= (2345) 


1(5) 


= (7234) 


1(4) 


= (7823) 


1(3) 


= (7812) 


1(2) 


= (4781) 


Id) 


= (4578) 



This data is generated by the function necklaceR according to: 



in[i] := necklaceR [{3,7,6,10,9,8,13, 12}] 
0ut[i]:= {{4, 5, 7, 8}, {1,4, 7, 8}, {1,2, 7, 8}, {2, 3, 7, 8}, 
{2, 3, 4, 7}, {2, 3, 4, 5}, {2, 4, 5, 6}, {4, 5, 6, 7}} 



(3.2) 



• nonSingularQ[permutation_,twistorDimension_:4 ] : tests whether or not a config- 

uration labeled by permutation has non- vanishing support for generic kinematical 
data — nonSingularQ returns True iff intersectionNumber [permutation}> 0. 

• parityConj ugat e[pcrmutation_ ] : if permutation labels a configuration C G G(k, n), 

then parityConjugate returns the permutation labeling the geometrically dual 
configuration C 1 - &G(n-k,n). For example, 

in[i] := parityConjugate [{3,5,6,7,8, 10}] 

parityConj ugat e@°/ 
out[i]:= {4,5,7,6,8,9} 

{3,5,6,7,8,10} 
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• permToGeometry[permutation_,removeVanishingQ_:Tr\ie ]: returns a (formatted) 
table of planes of various ranks spanned by consecutive chains of column-vectors 
of the configuration labeled by permutation; distinguished planes are highlighted 
in blue, and (for visual clarity) the option removeVanishingQ (True by default) 
causes any vanishing columns of the configuration to be suppressed. 

For example, consider the positroid configuration labeled by the permutation 

{3,7,6,10,9,8,13,12}, 

consec. chains of columns span 

(1) (2) (3) (4) (5) (6) (7) (8) P° 
(123) (678) P 1 ( SJ ) 

(56781) (81234) (3456) P 2 



For this, the function permToGeometry would produce: 

in[i] := permToGeometry [{3,7,6,10,9,8,13, 12}] 

Gutm== / (1) (2) (3) (4) (5) (6) (7) (8) \ 
(123)(34)(45) (56) (678) (81) 
\ (3456) (56781) (81234) / 

(Here, all the 'distinguished' maximal planes (a+1, . . . , cr(a)-l) have been high- 
lighted in blue — as the rest of the table follows form knowledge of these planes.) 

• permut at i onK[permutation_ ]: returns the "fc" associated with permutation. More 
explicitly, if C = (ci, . . . , c n ) is the positroid configuration labeled by permuta- 
tion, then permutationK returns rank{ci, . . . , c n }; as such, C G G(k, n). E.g., 

in[i] := permutationK [{3,7,6,10,9,8,13, 12}] 

Dut[l]:= 4 



• pref err edGauge[permutatioii_]: returns the lexicographically first non- vanishing 
minor of the configuration C G G(k, n); equivalently, if C = (c\, . . . , c n ) is the con- 
figuration labeled by permutation, then preferredGauge returns the lexicographically- 
minimal set of column-labels (oi, . . . , a k ) such that rank{c ai , . . . ,c ak } = k. 
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• randomCell[n_,ic_,d_,excIusionsQ_:True ]: gives the permutation label of a randomly- 

generated, d-dimensional positroid cell C G G(k,n). If the optional argument 
exclusionsQ is True (its default value), then randomCell tries to find a configu- 
ration for which all columns are non- vanishing (when possible). 

• rota.te[rotation_~\[permutation_]: returns the permutation labeling a positroid 

cell whose columns are rotated (positively) relative to that of permutation by ro- 
tation; more specifically, given a positroid cell C = (c\, . . . , c n ) G G(k, n) labeled 
by permutation, rotate returns the label of the configuration C = (cy, . . . , <v) 
where a' = a+rotation. 

• storeBoundaries: a global variable which by default is set to False; when True, 

boundary information is stored in memory so that boundary computations need 
not be repeated. 
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3.2 Positroid Coordinates and Matrix Representatives 



bridgeToMinors[pciTnutation_] : using the BCFW-bridge construction of coor- 
dinates for the configuration C labeled by permutation, bridge ToMinors ex- 
presses each BCFW-bridge coordinate aij directly in terms of the minors of C. 
For example, the {canonical) BCFW-bridge decomposition of the permutation 
{3, 7, 6, 10, 9, 8, 13, 12} generates the matrix representative: 



in[i] := permToMatrix [{3,7,6,10,9,8,13, 12}] //nice 

0ut[i]:= /lag -a$ -anas \ 
1 a 8 a 7 
1 a^+aQ 0304 —a\ 
\0 1 a 4 «2 / 

It is a highly non-trivial fact that such bridge coordinates can be expressed 
as ratios of monomials involving the minors of C. The explicit correspondence 
is given by the function bridge ToMinors, as in the following example: 



bridgeToMinors [{3,7,6,10,9,8,13, 12}] //r 

(1258) (1247) 
"l (1245) "2 -> (Y^) «3 - 



In[l] : = 
0ut[l] : = 



/nice 

(1247) (1245)(1278)(4567) 

,i^u, 2 (1245) 3 (1246)(1247)(4578) 

(1246)(1267)(4578) (1246)(4578) (1456)(2578) 

0:4 ~~ y 1245 (1278)(4567) 015 ~^ 1278)(1456) a% ~~ ^ 
(1478) , (1345) 

"8 ~^ (1245) " 9 



a 7 — > 



(1278) 



1246)(4578) 
2456) 



(14 5 6) 



• matrixCharts[permutation_] : gives a list of (distinct) matrix-representatives of 

the positroid configuration C labeled by permutation, using each of the (cyclically- 
distinct) columns as the 'minimal' column used in the construction of BCFW- 
bridge coordinates for C. 

• matrixToPerm[matrix_ ] : returns the permutation which labels the positroid cell 

represented by matrix. If matrix is given in terms of unspecified variables, then 
matrixToPerm assumes that all such take on generic values. 

• permToMatr ix[permutation_,transpositionSeheme_:0 ]: returns a matrix- representative 

of the positroid cell labeled by permutation given in terms of BCFW-bridge co- 
ordinates obtained using one of several possible bridge-decomposition schema — 
with the default scheme being "0", that corresponding to the canonical or 'lexi- 
cographic' scheme described in reference [1]. The possible decomposition schema 
are those described for transpositionChain. 
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' positiveQ[matn'x_] : returns True if the all ordered minors of the matrix are 
strictly non-negative; if matrix is parameterized by unspecified variables, then 
positiveQ returns True if matrix is positive (in the previous sense) when all its 
unfixed variables are given random, positive values. 

■ transpositionChain[permutation_,sciieme_:0 ] : returns a list {transpositionList , 
permutationList , seedGauge} containing the complete BCFW-bridge decom- 
position of permutation into 'adjacent' transpositions according to the scheme 
denoted scheme. (This data is obviously redundant, but such redundancy proves 
somewhat useful to have at hand.) The possible schema include: 



scheme 




(default) 


the canonical or 'lexicographic' decomposition scheme 


' 'cyclic' ' 


a scheme which attempts to decompose permutation into a sequence 
of bridges in a way which preserves cyclic symmetry (if any exists) 


1 


a scheme which correctly orients all (2n-4)-dimensional cells in the 
momentum-space Grassmannian 


-1 


a scheme which correctly orients all 4fc-dimensional cells in the 
momentum-twistor Grassmannian 



To illustrate the different transposition schemes, the lexicographic scheme 
( "scheme=0" ) would decompose the permutation {4, 7, 6, 9, 8, 11, 10, 13}, 

"0" (Lexicographic) Bridge Decomposition Scheme 



(12) 
(2 3) 
(34) 
(2 3) 
(12) 
(3 5) 
(2 3) 
(5 6) 
(3 5) 
(5 7) 
(2 5) 
(5 8) 



2 

i_ 

7 

4 

6 

6 

9 

7 

7 



8 
8 

9 10 
9 10 



3 

±_ 

6 

6 

4 

9 

6 

6 

8 

7 

7 
11 
11 
11 
11 



4 

9 
9 
9 



5 



6 
6 
11 

7 
10 

8 

13 



6 7 8 

\~ \~ 4" 

11 10 13 

11 10 13 

11 10 13 

11 10 13 

11 10 13 

11 10 13 

11 10 13 

11 10 13 

10 13 

10 13 

7 13 

7 13 

6 7 8 



BCFW shift 

c 2 c 2 + a 12 ci 
c 3 !->> c 3 + a u c 2 

C 4 H- c 4 + «i C 3 

c 3 ^ c 3 + a 9 c 2 
c 2 h-> c 2 + a 8 ci 
c 5 H- c 5 + a 7 c 3 
c 3 h-> c 3 + a 6 c 2 
c 6 >-> c 6 + a 5 c 5 
c 5 i-> c 5 + a 4 c 3 
c 7 h-). c 7 + a 3 c 5 
C5 ^ c 5 - a 2 c 2 

C 8 h-» C 8 + «1 c 5 
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This decomposition would give rise to the following representative of the cell: 
/l (08 + 012) (ag + aii)ag 0309010 ^ 

1 (a6 + a9 + ail) (06 + Og)oio (o 6 07-02) -0205 -0203 

1 oio (04 + 07) 0405 

yo 1 o 5 o 3 -oi/ 



However, the "cyclic" scheme would decompose the permutation according to: 
"cyclic" Bridge Decomposition Scheme 



2 

I 



3 
I 



4 

I 



5 
I 



6 

I 



7 

I 



7 
4 
4 
4 
4 
9 
9 
9 
9 
7 
7 
7 
2 



6 
6 
9 
9 
9 
4 
6 
6 
6 
6 
11 
11 
11 



8 11 10 13 

8 11 10 13 

8 11 10 13 

11 8 10 13 

11 8 13 10 

11 8 13 10 

11 8 13 10 

11 13 8 10 

11 13 10 

11 13 10 

6 13 10 
10 

13 6 15 



S 



BCFW shift 



c 2 h+ c 2 + a 12 Ci 

C4 !->• C4 + CK11C3 

ce ^ c 6 + «ioc 5 
c 8 c 8 + a 9 c 7 
c 3 H> c 3 + o 8 c 2 
C4 1 — y C4 + «7 C3 

C 7 H+ C 7 + 6 C 6 
C 8 H-). C 8 + 5 C 7 

c 2 1-+ c 2 + o 4 Ci 

C 5 ^ C 5 + « 3 C 3 

ce H- c 6 + a 2 c 5 

C2 !"+ C 2 — «l C 7 




which would result in the following matrix-representative: 



(l (04 + 012) 




\0 oi 



0408 04O8O11 
1 (07+011] 


Oi08 O1O8O11 



\ 

— 02 —0206 — 020609 

1 (a 3 + oio) a 3 06 a3a 6 a 9 

1 (a 5 + o 9 )y 



(3.5) 
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3.3 Drawing On-Shell (Plabic) Graphs and Left-Right Paths 



• plabicGr ajphlpermutation_,optionsList :def aultOpt ions] : draws a reduced plabic 

graph whose left-right path would be given by permutation. (Here, 'plabic' is 
used somewhat loosely, as the default behavior of plabicGraph is to draw graphs 
with monovalent and bivalent vertices when the inclusion of such is warranted.) 

There are many possible options for plabicGraph; the principle of these are: 



option value description 



chainOption * 





uses the lexicographic bridge decomposition scheme 


cyclic 


uses the bridge decomposition scheme "cyclic" 


±1 


uses the bridge decomposition scheme set by ±1 


rotation * 





constructs the graph using whatever bridge decomposi- 
tion scheme is specified, but where particle '1' is consid- 
ered minimal 


a 


same as above, but with cyclic ordering beginning with 
particle (1— a) 


LRpaths * 


{} 


does not show any left-right paths 


{a, . . . ,bj 


draws the left-right paths which start at legs a, . . . ,b 


"All" 


draws all left-right paths 


directed * 


False 


results m an undirected graph 


True 


results in a directed graph whose perfect orientation fol- 
lows from the BCFW bridge-decomposition scheme used 


edgeLabels * 


False 


does not label any edges of the graph 


True 


laueis uie cur vv-DiiQge euges uy tne Diiuge cooiQinates 
associated with each 


faceLabels * 


False 


does not label any faces of the graph 


True 


labels each face of the graph with a label 


"A" 


labels the faces of the graph with A-variable labels 


<< x> , 


labels the faces of the graph with X-variable labels 


showRemovable * 


True 


shows monovalent vertices attached to legs which are 
self-identified under permutation 


False 


removes boundary legs which are self-identified under 

permutation 


orientation * 


1 


draws the external legs with clockwise ordering 


-1 


draws external legs with counterclockwise ordering 


bipartiteQ * 


False 


allows for trees of same-colored vertices 


True 


collapses all same-colored trees giving a bipartite graph 
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Here, each default option has been marked with a In addition to these op- 
tions, there are many which control stylistic and aesthetic details of the graphs 
generated by plabicGraph; the principle among these include: 



option 




value 


description 






n 


(~\y QXTlTC TrlO (TV Q T~\ rl \~XT1 T rl T~\ Q T"T 1 C 1 T T CVtl7Q Tn f (ID TT1 1 CI CI \ d 1 DTI 

LlidWo Lllc ^IdUll W1L11 pdi LlLlc 1 LUWdiU. Lllc 1111 LlU.lt; - It; 1 L 








vot n t pc flip orann \~\~\r r) ra H 1 53 n c: in tnp o I no la no qp Hi TPftinn 

iULaLCO Lilt- til aUll Uy U ICLUldllo 111 Lilt; ULULsFl W IjO l> Ull CL LlUll 






fill uu 


IdUclo lilt; tj7vLt;i lldl ItJtia 1 1 } . . . 3 11 r 






fa bl 


la nplc; Trip P'VTPTnal 1po*c; J a \\ I 

IdUClO ullC CALCI 11CL1 1C&Q l CL j . . . 3 U f 






fiR 


Filfippc; fnp P'VTPTTial FifiTfiplp lanplc; nf a nistpTifP ot fiR 

UldL-Cu L11C- LC<1 11CL1 IJCll LHjlC IdUClO CL u CI U.10 uCLllLC Ul \J • w w 

from the boundary 
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vertexSize 




0.325 


draws all vertices with size of . 325 units 


lineThickness 




3.5 


draws edges with AbsoluteThickness [3 . 5] 


LRpathTh i ckne s s 




4 


draws left-right paths with AbsoluteThickness [4] 


LRpathDistance 




0.275 


draws left-right paths with a distance of . 275 units 
from the graph's edges 


LRarrowHeadSize 




0.0655 


sets the size of left-right path Arrowheads to be . 0655 


edgeArrowSize 




0.05 


sets the size of directed-edge Arrowheads to be . 05 


outerCircle 




True 


draws a disc at the boundary, enclosing the graph 



Wherever sufficiently obvious, we have left as implicit the possible alternative 
settings allowed these options. 



For a complete list of options for plabicGraph — together with the default 
value for each — consult the global variable "defaultPlabicGraphOptions" ; users 
should find most detailed features of the output to be pliable through simple 

experimentation with the options. 

• plabicGraphData[permutation_, transpositionScheme_:0 ]: given the permutation 

label for a positroid configuration a some transpositionScheme, plabicGraphData 
returns the pair {edgeList ,f aceList}, where edgeList lists each (directed) 
edge i— >j (carrying a weight w) as a triple {i,j,w}; and f aceList lists the 

vertices along each face of the graph, listing them with clockwise ordering. 

• resetGraphDef aults:a function which automatically resets all the default graph 

options for the function plabicGraph. (By default, plabicGraph remembers any 
options explicitly set until the length of permutations being drawn is changed.) 
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3.4 Physical Operations, Kinematics, and Scattering Amplitudes 



• bcfwPartitions[ri_,.k_]: returns alist of pairs {(til, (tir, /cr)} which are bridged- 

together to compute the n-point N^^MHV tree-amplitude according to 

4 k) = £ A t ] (8) A nn- (3- 6 ) 

(n L ,k L ) BCFW 

• bcf wTermNames[i2_,I_] : returns a formatted list of terms appearing in the BCFW 

tree-amplitude where MHV and MHV amplitudes have been collected together. 

in[i]:= bcf wTermNames [6 , 3] 

bcfwTermNames[14,7] [[9598]] 

out [i] : = {4 3) ® 4 13 , 4 2) ® 4 2) , 4 2) ® 4 2) } 

( (4 2) ® (4 2) ®4 2) ) ) ®4 } ) ® (4 2) ® ( (4 2) ® 4 2) ) ® 4 X) ) ) 



• generalTreeContour[a_:0, b_:0,bridgeChoice_:Q ] ln_,k_ ] : returns a list of permu- 

tation labels for the positroid cells which together give the n-particle N( k - 2 )MHV 
tree-amplitude obtained using the white-to-black BCFW-bridge attached to legs 
(In) for bridgeChoice '0' (the default) or (nl) for bridgeChoice '1', and for 
which the lower-point amplitudes appearing in the recursion have been recursed 
(using the same bridge-choice) attached to legs (1+a, n^ + a) and (1 + b, riR+b) 
of the left and right amplitudes, respectively. 

• identitySigns[permList_] : returns a list of ±1 of the relative signs needed for an 

identity involving the cells labeled by the permutations in permList. E.g., 

in[i]:= identitySigns[{{3,2,4,5,6,7},{2,4,3,5,6,7},{2,3,5,4,6,7}, 
{2,3,4,6,5, 7} ,{2, 3,4,5, 7, 6} ,{1,3,4,5,6,8}}] 

0ut[l]:= {1,-1,1,-1,1,-1} 



• nRatioContour[n_,ic_]: produces the same output as: 

permToResidue/ OdualGrassmannian/ @treeContour [n,i] . 

• nTreeContour[n_,]f_] : produces the same output as: 

permToResidue/ OtreeContour [n,.k] . 
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' permToResidue[permutation_ ] : for a permutation labeling either a (2n-4)-dimensional 
cell C in the momentum-space Grassmannian or a 4fc-dimensional cell in the 
momentum-twistor Grassmannian, permToResidue uses the globally defined mo- 
mentum twistors Zs and the corresponding (or alternatively-defined) global, 
spinor variables Ls and Lbs (A and A, respectively), to find the isolated point (s) 
C*6C which solve the kinematical constraints and returns a pair (or list of pairs 
if more than one solution exists) where /x* is the positroid measure 

evaluated at the point C* which solves the kinematical constraints. 

The on-shell function corresponding to a positroid cell labeled by a, when 
evaluated at whatever kinematical data is given, would be given by, 

f a = lSx6 kxA (C*.rj), or f a = t i*x5 kx4 (C*- V ), (3.7) 

if a labels a cell in the momentum-space Grassmannian or the momentum-twistor 
Grassmannian, respectively. 

■ positiveZs[nPartic!es_] :it is sometimes convenient to evaluate analytic expres- 
sions involving spinor-helicity variables or momentum-twistors using explicit 
kinematical data; under such circumstances, there are some conveniences af- 
forded by using "well-chosen" kinematical data. 

Reasons for preferring one choice over another include: having all Lorentz 
invariants be integer-valued and relatively small; having all dual-conformal cross 
ratios positive (so as to avoid branch-ambiguities when evaluating the polyloga- 
rithms that arise in scattering amplitudes at loop-level); and possibly to have all 
Lorentz-invariants be distinct (either to help reconstruct an analytic expression 
or to avoid 'accidental' cancelations). Of these, the following momentum-twistors 
meet the first two desires spectacularly: 

/l 1 1 1 ■■■ Co) \ 

2 3 4 5 ■■■ 

3 6 10 15 • • • ( n + 2 ) 
^4 10 20 35 • • • ( n + 3 ) J 



Zs 



(3i 



The function positiveZs [16] is evaluated when the positroids package is 
first loaded, allowing amplitudes involving as many as 16 particles to be evaluated 
without specific initialization. 

randomTreeContour[ri_,i_] : returns the list of permutation labels for cells occur- 
ring in the BCFW tree-amplitude formula, where the lower-point amplitudes 
have been recursed using randomly-chosen legs. 
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• setupUsingSpinors[iambdaList_,iambdai3arList_]:sets up the global variables 

Ls and Lbs for A and A, respectively, and defines the global (nx4) matrix Zs for 
momentum twistors for use in numerical evaluation. 

• setupUsingTwistors[twistorList_]:sets up the global (nx4) matrix Zs encoding 

the momentum-twistor kinematical data, and defines the auxiliary variables Ls 
and Lbs for A and A, respectively. 

• superComponent[component ] [superFunction_ ] : for the purposes of the positroids 

package, a superFunction must given by a pair {/, C} — an ordinary function 
f(l, . . . , n) of the kinematical variables times a fermionic 5-function of the form, 

S M (C-rf)=f[l (ay--a k )vl r -vi\, (3.9) 

1=1 Loi<— <OJ J 

where C is an (n x fc)-matrix of ordinary functions, and for each a — 1, . . . , n, T) a 
is a fermionic (anti-commuting) variable. To be clear, we consider each particle 
Grassmann coherent state of the form, 

I \ I \ ' — J I \ 1 I \ 1 -^r, — 1 i \ : — 1 * — 1' — '\ — 4 i \ 

| a) = \a} {} + r] a \a) {I} + -r] a r] a \a) {IJ} + ^r] a r] a r] a \a) {I JK} + VaV a V a Va W{i,2,3,4} 5 
and if we use r a to denote the i?-charge of the a th particle according to, 



field 


helicity 


/2-charge (r a ) 


short-hand for r a 




+1 


{} 


P 


\a){i} 


+ \ 


{1} 


P/2(« {4}) 


\ a ){i,J} 





{I>J} 




\ a ){I,J,K} 


l 

2 


{I,J,K} 


m/2(^ {1,2,3}) 


l a ){l,2,3,4} 


-1 


{1,2,3,4} 


m 



then superComponent[ri,. . . ,r„] [superFunction] returns the component function 
of superFunction proportional to, 

niK- (3 - io) 

a=l Ier a 

— that is, the component-function involving the states: 

|l) ri ••»,„■ (3.11) 

For example, the "(— , +, — , +, — , +, — , +)" component of the 8-particle N 2 MHV 
amplitude proportional to, 

would be extracted be obtained by evaluating, (compare with e.g. [23,24]): 

in[i] := Total [superComponent [m,p,m,p,m,p,m,p] /@ 

permToResidue/ OtreeContour [8,4]] 
out[i]:= 908 416 
~ 39 375 
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• terms InBCFW[n_,.k_/_:0 ]: gives the number of (non-vanishing) terms generated 
by the BCFW-recursion for the (?-loop, n-point N^ k_2 ^MHV amplitude. (£ must 
be either or 1 as the number of terms is scheme-dependent beyond 1-loop). 

in[i]:= termsInBCFW[6,3] 

termsInBCFW [6,3,1] 

Out[l]:= 3 
16 



• treeContour[n_,i_ ] : returns the list of permutation labels for positroid cells which 
together give the n-particle N^ k_2 ^MHV tree-amplitude, using the default recur- 
sion scheme; treeContour [n,k~\ is equivalent to generalTreeContour [0,0,0] ln,k] . 

in[i]:= treeContour [6 , 3] 

outci]:- {{4,5,6,8, 7, 9}, {3, 5, 6, 7, 8, 10}, {4, 6, 5, 7, 8, 9}} 
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3.5 Aesthetic and General Purpose Functions 

• explicif y[expression_,positiveQ_.True ] : picks random (integers) for all variables 

occurring in an expression; if positiveQ is True (its default value), then the 
random assignments are taken to be positive. (If expression includes angle- 
brackets (•••), explicify will evaluate these assuming random kinematical data.) 

• exportToPDF[ii]eiVame_ ] \_figure_ ] : saves a PDF version of figure to the file file- 

name 4 using outlined fonts (and with other minor processing). 

• mod2[objectList_ ] : returns the elements of objectList which occur an odd number 

of times in objectList. This can be useful, for example, if one wants to explicitly 
verify that d 2 = mod 2: 

in[i] := mod2 [Join@@(boundary/@boundary [randomCell [8,4, 12]] )] 

Out[l]:= {} 



• nice[expression_] : formats expression to display 'nicely' by making replacements 

such as ab [x- • • y] \-> (x---y), a[l]i— >■ a±, etc., and by writing any level-zero 
matrices in MatrixForm. 

• nice! ime[timeInSeconds_~\: converts a time measured in seconds timelnSeconds, 

to human-readable form. For example, 

in[i]:= niceTime[299 792 458] 

niceTime [3 . 1415926535] 
out[i]:= 9 years, 182 days 

3 seconds, 141 ms 

• random[objectList_ ] : returns a random element from (the first level of) objectList. 

• randomSubset[subsetLength_] [objectList_] : returns a randomly-chosen subset of 

length subsetLength from among the list objectList. 

• timed[expression_] : evaluates expression and prints a message regarding the time 

required for evaluation. 



4 The file is saved to the same directory as Export [] : either NotebookDirectory [] or Directory [] 
(if the former does not exist). 
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